Magnetic traveling-stripe- forcing: 
enhanced transport in the advent of the Rosensweig instabiUty 
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A new kind of contactless pumping mechanism is realized in a layer of ferrofluid via a spatio- 
temporally modulated magnetic field. The resulting pressure gradient leads to a liquid ramp, which 
is measured by means of X-rays. The transport mechanism works best if a resonance of the sur- 
face waves with the driving is achieved. The behavior can be understood semi-quantitatively by 
considering the magnetically influenced dispersion relation of the fluid. 
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I. INTRODUCTION 



Ferrofluids are colloidal dispersions of magnetic 
nanoparticles, which show super-paramagnetic behavior 
This offers the advantage to control and facilitate 
their flow by time-dependent magnetic fields. For ex- 
ample in alternating fields a viscosity reduction was pre- 
dicted for Poiseuille flow Q and thereafter experimentally 
confirmed 0, More spectacularly, in rotating fields a 
spontaneous spin up of surface fiow was observed in a 
beaker [51] , and quantitatively analyzed in a circular duct 



An intriguing alternative to rotating fields is to apply 
a spatio-temporally modulated field. Here one has to 
discern two transport mechanisms. One is "ferrohydro- 
dynamic pumping" which was lately studied for a closed 
channel geometry [7, 8]. It relies on the phase shift be- 
tween the external field and the magnetization and works 
best for frequencies in the range of the inverse Brown- 
ian relaxation time, i.e. in the kHz-regime. The second 
mechanism relies on traveling surface waves, driven by a 
traveling magnetic field. This was realized by an array of 
solenoids acting on a free surface of ferrofluid [9*1, and dis- 
cussed in connection with the effect of peristaltic pump- 
ing As we will show in this article the driving caused 
by a free surface gains momentum if the traveling field 
is superimposed by a constant field oriented normally to 
the fluid layer. 

In such a vertically oriented magnetic field, a plain 
layer of ferrofluid disintegrates into a hexagonal array of 
crests if a critical threshold He is overcome - the so called 
Rosensweig instability [11]. Approaching He from below, 
a vanishing of the phase velocity Vp of surface waves was 
observed, and explained by the non-monotonic dispersion 
relation [12]. In this regime, a traveling-stripe-forcing of 
the magnetic induction was predicted to excite a reso- 
nance of the wave amplitudes if the driving velocity co- 
incides with Vp [l3| . This was recently verified by exper- 
iments [Hi. As we show in the following, this resonance 
is correlated with a maximal net force - as indicated by 
a liquid ramp - on the magnetic fluid. 
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FIG. 1: Sketch of the experimental setup 




FIG. 2: The crosses indicate the magnetic field as a function 
of time. The horizontal gray line corresponds to the applied 
field of the Helmholtz coils fioHo as explained in the text. The 
black line represents a fitted harmonic function. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

Our setup is sketched in Fig.[T] A Helmholtz pair of 
coils provides an applied field Hq along the vertical di- 
rection by means of a stabilized electrical current. In the 
center of the coils a container machined from Perspex^*^ is 
located. It contains a box shaped cavity with a length of 
120 mm, a width of 100 mm and a height of 25 mm. This 
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FIG. 3: Modulation amplitude vs. the induction impressed by 
the coils. The experimental data (black crosses) are recorded 
for a distance of z = 2 mm between the Hall probe and the 
ground of the container. Gray line: linear fit. 

cavity is filled with 30 ml of ferrofluid EMGgog (Tab.[l|. 
About 5 mm below the bottom of the cavity, a motor 
driven conveyor belt is located. It harbors a grid of rods 
made from welding wire (0 = 2 mm) with a spacing of 
Ag = (9.3± l)mm. 

Radiation emitted by an X-ray tube mounted 1.5 m 
above the container is permeating the ferrofluidic layer 
and recorded by an X-ray detector (16 bit). The absorp- 
tion pictures serve to reconstruct the surface topography 
of the liquid layer. The measurement of the heig ht is 
calibrated by means of a ferrofluidic ramp [l^ Il6l |. To 
eliminate the shadows cast by the rods we measure the 
time averaged absorption of X-rays with the moving belt 
included. 

The belt leads to a spatial modulation of the magnetic 
field. The modulation of this field was measured within 
the empty cavity at one fixed position 2 mm above the 
ground of the container with the belt in motion. A result 
for a belt velocity vq of 0.75 ^ and an applied field of 
IIqHq = 22.29 mT, which has been measured by the Hall 
probe in the absence of the grid, is depicted in Fig. [51 
The fit to a harmonic function demonstrates that the 
magnetic field can be modeled as a harmonic wave 

Hg ^ HGfi + AHg sin {u}t~ kax). (1) 

Here i?G,o denotes the mean value and AHg,o the ampli- 
tude of the modulated field. The frequency w — 27rv^|14j 
is controlled by the motor (Mattke MDR-230). The 
wavenumber kG — is set by the grid spacing, and 
has been chosen to be close to the critical wavenumber 
of the Rosensweig instability. The deviation between the 
value of the applied field Hq as indicated by the gray line 
and the mean field Hg.o is explained by a focusing of the 
magnetic field due to the grid. 

Fig. [3] shows the modulation amplitude AHg as a func- 
tion of the applied magnetic field Hq. The fitted gray 
line indicates a linear dependence, thus demonstrating 
that the soft ferromagnetic material chosen for the rods 
is well suited for the experiments. 
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FIG. 4: AHg measured for different distances z at hoHq = 
22.29 mT (black crosses), gray line: exponential fit. 

The modulation AHg decreases with increasing dis- 
tance from the rods. This is demonstrated experimen- 
tally in Fig.m The gray line is the fit to 

AHG{z) = AHG{0)e-'^^^, (2) 

because we expect an exponential decay with the decay 
length given by the grid spacing when assuming a two 
dimensional field. The fit yields //q Ai?G(0) = 0.62 mT 
for the modulation at the bottom of the cavity. 

The ferrofluid EMGgog is made of magnetite particles 
dispersed in kerosene. Its nonlinear magnetization curve 
is plotted in Fig. [51 We measured M{H) using a flux- 
metric magnetometer consisting of a Helmholtz pair of 
sensing coils with 6800 windings and a commercial inte- 
grator (Lakeshore Fluxmeter 480). The sample is held in 
a spherical cavity with a diameter of 12.4 mm in order to 
provide a homogeneous magnetic field inside the sample 
with a demagnetization factor of ^ . The gray line is a fit 
using the Langevin function [l| 

M ^ Ms (^coth {/3H)~^y (3) 

This function was derived for dilute monodisperse col- 
loidal suspensions and shows visible deviations from the 
data. In contrast, the black line displays a fit of a model 
for dense polydisperse magnetic fiuids, namely Eq. (32) 
of Ref. [17]. This model assumes a F- distribution for the 
particle diameter d 

where do ~ 1.4 nm and a — 3.8 are the obtained fit 
parameters, corresponding to a mean particle diameter 
d — 6.6 nm. The fit also provides the volume fraction of 
magnetite to be 4.0%. 

Furthermore we characterize our ferrofiuid by subject- 
ing it to a homogeneous magnetic field, i.e. we omit for 
this measurement the driving belt. The plain layer be- 
comes unstable at a critical field of /io^c.meas = 21 mT. 




FIG. 5: Nonlinear magnetization curve M{H). Black crosses: 
measured data, black line: fit using the Ivanov model [l7l |. 
gray line: fitted Langevin function. 



Like in Refs. 

this value was obtained by fitting the 
measured spike amplitudes with a scaling law provided 
by theory [1^. 

Figure [n](a) shows the surface topography of the mag- 
netic liquid as reconstructed from the X-ray absorption, 
for a driving velocity of = 4.1 ™ to the l.h.s. We 
averaged over 1300 subsequent images with an exposure 
time of 0.4 s each. Obviously the liquid builds up a ramp 
in the container. A longitudinal cut of the ramp is pre- 
sented in Fig.[S](b) by black crosses. Each data point 
stems from an average along the y-dimension selecting 
the innermost 36.4 mm of the vessel. The gray symbols 
give the longitudinal cut for a driving to the r.h.s. The 
liquid ramp can be approximated by 



h{x) = tan {j)x + ho- 



(5) 



as indicated by the fitted lines. The difference of the 
inclinations is measured by the angle 27. This value is 
representative for the pressure gradient created by the 
magnetic pump. 



III. AN ESTIMATE OF THE PRESSURE 
GRADIENT 

The build up of the liquid ramp is caused by the Kelvin 
force, since we apply an inhomogeneous magnetic field to 
magnetized matter. The system can be described by the 
ferrohydrodynamic Bernoulli-equation fll| 



dv 



= -Vp + ^ioMVH + -qV'^v + pg. (6) 



As a model assumption, we assume our system to contain 
no fiow, thus all the terms with v vanish. Further, we use 
a one-dimensional ansatz, and get 



dp OHg 

— — tlQlUcS^ ■ 

ox ox 



(7) 



The strength Hq of the external magnetic field in har- 
monic approximation is given by Eq. ([T|). The effective 
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FIG. 6: A liquid ramp: (a) Three dimensional reconstruction. 
The height difference between three subsequent black lines is 
1.0 mm. (b) Averaged height profiles for forcing to the l.h.s 
(black) and to the r.h.s. (gray). The solid lines represent 
fits by Eq. ([5]) within the range between x = 60 mm and x = 
120 mm. 



magnetization is defined as 

Meff (a;) = M 



h{x) 



(8) 



where ho denotes the average height of the ferrofiuid layer, 
and M is assumed to be constant. Assuming a harmonic 
wave for the height modulation 



h{x) = /iQ + A/i sin {ut — kcx + (/)), 



(9) 



where (/) denotes the phase difference with respect to the 
magnetic field, one gets 

Mcff{x) = Mcffap sin {ut - kax) 

- Mcff,opcos(a;t - fcca;) + Mcff,o- (10) 

Here A/off . ip denotes the part of the magnetization which 
is in phase with the external field according to Eq. ([1]). In 
contrast, Moff_op is out of phase with Hq, but in phase 
with the gradient 



dHGix,t) 
dx 



= —kAHc cos {bjt — kox) 



(11) 



of the external field. Inserting Mcs from Eq. (jlOp into 
Eq. ([7]) and averaging over one period of the external field 
leads to a mean pressure gradient 



dp 
dx 



/oVoAfcff'^dx 



(12) 
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along the x direction. This pressure gradient is experi- 
mentally determined with the observed inclination 7 of 
the surface: 



dp 
dx 



P5tan7. 



(13) 



IV. EXPERIMENTAL RESULTS 

In order to characterize the pressure difference created 
by the pump, we investigate its dependence on the exter- 
nally applied field Hq and the belt velocity vq- Figure [7| 
presents the results for ten different values of Hq. For 
each field strength, vq was raised from to 11.2 ^ in 30 
steps. 

Although the width and the shape of the 7(z;)-curves 
changes with the applied field, in each case the pressure 
difference exhibits a maximum at a certain driving veloc- 
ity i'm(-ffo), e.g. for /j,oi?o = 18.95 mT 7 reaches a max- 
imum at 7.29^. With increasing iJo this maximum is 
shifted to lower velocities. In order to determine the po- 
sition of the maxima (wm), we calculate the barycenter of 
those data being larger than 75% of the maximal value 
of 7. 

In Fig.[5]a) the outcome is marked by black crosses. 
The gray curve denotes a model presented in the next 
section. The maximal surface inclination is extracted by 
fitting parabolas to the upper part of the curves. The 
result is plotted in Fig.[5]b). It demonstrates a monotonic 
increase of the efficiency of the pump with the applied 
magnetic field, which seems to saturate at about 30 mT. 
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FIG. 7: Selection of surface inclinations versus the driving 
frequency for different fioHo (denoted by the right ordinate 
in mT), the curves are shifted equidistantly upwards to avoid 
intersections. The data points are connected by splines to 
guide the eye. 



ity 'At k = ka one yields 



V. MODELING THE RESONANCE 

The resonance of 7 for a specific Vm, as depicted 
in Fig. [3 resembles the resonance of the amplitudes 
of surface waves below the onset of the Rosensweig 
instability [H, [l3|- The resonant phase velocity Vp{Ho) 
can be captured by the dispersion relation for an inviscid 
magnetic layer of infinite depth fl^. 

According to Eq. (fT2| . we expect a monotonic relation 
between the wave amplitudes and the inclination 7. The 
function VmiHo) in Fig.[5]a) can be modeled in a certain 
regime by utilizing the dispersion relation for small sur- 
face wave amplitudes 



uj'^{k) ^ g^k + - 



Mo 



p r + 1 



(14) 



It can be deduced from Eq. (36) in Ref. [20'], as detailed 
in Ref. [21j . The effective permeability r is the geometric 
mean r = V(l + Xta)(l + Xch), where Xta = (dM/dH) 
denotes the tangential susceptibility, and Xch = (M/H) 
the chord susceptibility [1]. For the resonant phase veloc- 
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(15) 



Using the jump condition for the field at the bottom of 
the container i?cxt = + M{Hi) and assuming that 
i?oxt is approximately given by iJo, one can obtain Vp{Ho) 
numerically. The outcome is marked in Fig.|S]a) by the 
solid line. It has been calculated using the magnetization 
curve as presented in Fig.O It shows a good agreement 
with the measured data only for small values of the field, 
because it is based on the assumption of small amplitudes 
of the surface waves. The intersection of this line with 
the X-axis at /io^^c.caic = 21.35 mT denotes the estimate 
for the onset of the Rosensweig instability, where surface 
deformations spontaneously form without any external 
modulation. This critical value has been independently 
measured (see sec.|lT]and Tab.|I| and matches within the 
experimental uncertainty. 

The maximum of the pressure gradient has experimen- 
tally been observed for driving at Vm as shown in Fig.|8]b). 
The phase shift between Hq and M^s is assumed to be 
■| for this resonant driving, in analogy to a driven har- 
monic oscillator. Using Eqs. and this assump- 
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FIG. 8: a) Resonant phase velocity versus applied magnetic 
field. The black crosses give the experimental results. The 
gray line stems from Eq. (|15p ; b) Angles of maximal surface 
inclinations over the applied magnetic field. Black symbols 
are experimental data points. The black line just connects 
the points via a spline fit to guide the eyes. The gray circles 
are calculated according to Eq. ([16]). 



tion leads to 



pg \dx 



2pgho 



(16) 



We measured Ah with a laser reflection method, as 
described in Ref. [l3|. When using this method, we as- 
sume a sinusoidal modulation of the surface, it is thus 
restricted to small values of the magnetic field Hq . There- 



fore we are restricted to values well below TJc.noniin- The 
corresponding data points for four measured values of 
Ah are indicated in Fig.[5]b) by circles. The data only 
agree within about 50%. This is partly due to i) the 
flutter in Ac, ii) the geometrical uncertainty connected 
with the laser reflection method, iii) the variation in Hq 
along the ramp, iv) the variation of the magnetization 
and the field modulation within the layer thickness and 
v) the neglecting of the higher harmonics of h{x,t). All 
these uncertainties are taken into account by the error 
bars. Moreover the crudeness of the model should not be 
brushed under the carpet. In particular it does not take 
into account any backflow effects and assumes a strictly 
two dimensional geometry. 



TABLE I: Parameters measured for the ferrofiuid, EMGgog 
Lo t H030308A from Ferrotec Co. 

density p 0.9945 g/cm^ 

surface tension a 23.37 mN/m 

initial susceptibility xo 0.95 

viscosity -q 3.3 mPas 

critical wavelength \c = 2iy ./ -2- 9.89 mm 

° Y P90 

calculated critical field /loi^c.caic 21.35 mT 



measured critical field poH, 



c.mcas 



21 mT 



VI. CONCLUSION 

In this paper we have demonstrated a new kind of 
pump for magnetic fluids which is based on a resonance 
phenomenon in the advent of the Rosensweig instabil- 
ity. Our driving relies fundamentally on the excitation 
of surface waves. It is thus restricted to open channel 
geometries, and differs substantially from "ferrohydrody- 
namic pumping" as described in Ref. [8] which works in 
closed channels. That method, however, is most effective 
for frequencies comparable to the inverse Brownian re- 
laxation time (pa 10 kHz), while our driving mechanism 
favours low frequencies (w 10 Hz). 

We have characterized the pump under no-load con- 
ditions by measuring the static pressure gradient via a 
radioscopic method. Under variation of the driving ve- 
locity, the pressure gradient shows a resonance. The res- 
onant phase velocity decreases with the applied magnetic 
fleld. This can be understood by calculating the phase 
velocity of surface waves from the dispersion relation, tak- 
ing into account the nonlinear magnetization curve. The 
pressure gradient is observed to increase with the fleld, 
which is captured by a two dimensional model derived 
from the ferrohydrodynamic Bernoulli equation assuming 
a static equilibrium. This ansatz leads to a transport ef- 
fect, whose strength grows linearly with the surface mod- 
ulation, the peristaltic pumping on the other hand grows 
with its square. A refined model should include both 
effects and would be a task for numerical investigations. 

In future work, the spatio-temporal driving of the liq- 
uid could be achieved in different ways. While in the 
present study we used a moving belt, harboring iron rods, 
to drive the liquid, alternatively the driving may be real- 
ized by a network of current carrying wires without any 
moving parts. This would allow to control the modula- 
tion amplitude and the bias of the magnetic field inde- 
pendently. 

The pump has been tested under no-load conditions. 
The opposite approach would be to characterize the 
pump under maximum-load conditions. Therefore, the 
input and output of the pump must be connected by an 
open duct with minimal flow resistance, similar to the 
setup described in Ref. @. 

Note that our resonant driving may be exploited as 
well to drive non-magnetic liquids via the magnetically 
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driven peristaltic motion investigated in Ref. [2^. 

So far we have focused our measurements on the reso- 
nant behavior below the formation of surface undulations. 
It remains to be investigated how the spatio-temporal 
driving interacts with the secondary instabilities leading 
to more complicated surface patterns (23l |. 
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